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Abstract. The development of the ring current ions in the inner magnetosphere during the main 
phase of a magnetic storm is studied. The temporal and spatial evolution of the ion phase space 
densities in a dipole field are calculated using a three dimensional ring current model, considering 
charge exchange and Coulomb losses along drift paths. The simulation starts with a quiet time 
distribution. The model is tested by comparing calculated ion fluxes with Active Magnetospheric 
Particle Tracer Explorers/CCE measurement during the storm main phase on May 2, 1986. Most 
of the calculated omnidirectional fluxes are in good agreement with the data except on the dayside 
inner edge (L < 2.5) of the ring current, where the ion fluxes are underestimated. The model also 
reproduces the measured pitch angle distributions of ions with energies below 10 keV. At higher 
energy, an additional diffusion in pitch angle is necessary in order to fit the data. The role of the 
induced electric field on the ring current dynamics is also examined by simulating a series of 
substorm activities represented by stretching and collapsing the magnetic field lines. In response to 
the impulsively changing fields, the calculated ion energy content fluctuates about a mean value 
that grows steadily with the enhanced quiescent field. 


1. Introduction 

During the main phase of a magnetic storm, energetic 
plasmas are injected on the nightside from a boundary near the 
geosynchronous orbit. In response to the convection electric 
field, these particles drift inward and are trapped by the 
geomagnetic field and form the storm time ring current. The 
development of the ring current during the main phase of a 
magnetic storm has been studied. Wolf et al. [1982] applied 
the Rice convection model [Harel et al ., 1981] to study the 
early main phase of the storm of July 29, 1977. They found the 
magnetospheric convection, including effects of shielding, 
was sufficient to inject enough plasma sheet plasma deep into 
the magnetosphere to form a storm time ring current. Chen et 
al. [1994] used a guiding center simulation to model the main 
phase of storms as a sequence of substorm-associated 
enhancements in the convection electric field. They found that 
for storms with a main phase of about 3 hours the ring current 
enhancements are mainly associated with ions injected from 
open trajectories to closed ones. The ring current is augmented 
also by diffusive transport of higher-energy ions (E > 160 
keV) for storms having longer main phases (> 6 hours). They 
also found that transport alone would not account for the entire 
decrease in Dst typical for a major storm but an increased 
boundary value of the phase space density was necessary. 

Copyright 1996 by the American Geophysical Union. 

Paper number 96JA01274. 

0 1 48-0227/96/96J A-0 1 274$09 . 00 


In this paper, some problems of convection, collisional 
loss, and pitch angle diffusion effects on the ring current 
development during the storm main phase are addressed. The 
main phase of the magnetic storm on May 2, 1986, is 
simulated using a three-dimensional ring current model 
described by Fok et al. [1993, 1995a]. The model assumes a 
dipole magnetic field and the convection electric field of 
Volland-Stern, parameterized by the Kp index [Volland, 1973; 
Stern , 1975; Maynard and Chen , 1975], is employed. The 
special features of the model are the consideration of ions with 
arbitrary pitch angle distribution and the inclusion of 
Coulomb interactions with the thermal plasma in the 
plasmasphere, This model is developed not only for studying 
the ring current dynamics but it is also a useful tool to examine 
the interactions between the ring current and the background 
neutral and plasma populations [Fok et al ., 1993, 1995b]. The 
description of the model and the results in simulating the 
recovery phases of a moderate and a major storm are described 
briefly in the next section. In the following, only the 
evolution of H + ions during the development of the storm is 
discussed. The simulation starts with a quiet time distribution 
[Sheldon and Hamilton , 1993] at L shells between 2 and 6.5. 
The nightside energy spectra of H + fluxes measured by the 
Active Magnetospheric Particle Tracer Explorers 
(AMPTE)/CCE at L ~ 6.75 at every 16 hours serve as the model 
boundary conditions. The measured omnidirectional fluxes and 
pitch angle distributions (PADs) inside the injection boundary 
are compared with model predictions. 

The convection electric field model is a crucial factor in 
determining the particle drifts and thus the global ring current 


15,311 



15,312 


FOK ET AL.: RING CURRENT DEVELOPMENT DURING STORM MAIN PHASE 


dynamics. Usually, the convection field can be separated into 
an electrostatic component and a substorm-induced 
component. In this study, the different roles of the two 
components of the electric field on the ring current 
development will be examined. The ring current model will be 
run with an induced electric field in addition to the Volland- 
Stern-type convection, and the results will be compared with 
those from calculations in which the induced field is not 
included. The algorithm we use to construct the time-dependent 
induced field is different from those presented in previous 
works, and it will be described in section 6. In the following, 
average flux means the average differential intensity (in 
s'^cm^sr^keV' 1 ) over pitch angles, except where another 
definition is specifically given. 

2. The Model 

The 3-D ring current model solves the following bounce- 
averaged kinetic equation of the phase space density of the 
ring current ion species, considering charge exchange and 
Coulomb losses along the drift path [Fok et al, 1993, 1995a]: 

where f s (R 0 , <p , M, K) is the average phase space distribution 
function of species s along the field line between mirror 
points, R 0 is the radial distance at the equator, <p is the 
magnetic local time, M is the magnetic moment, v is velocity, 
G s is the cross section for charge exchange of species s with 
the neutral hydrogen, n H is the hydrogen density, and M is the 
rate of change of M due to Coulomb interactions with the 
thermal plasmas. K is defined as [Mcllwain, 1966], 

K = j/(Sm s M) m ( 2 ) 

and it is an invariant in the case of zero parallel electric field 
[. Roederer , 1970]. The notation <x> is the bounce-averaged 
value of x. Since we are considering particles with bounce 
periods much shorter than the decay lifetimes, f s is assumed to 
be constant along the field line, and thus f s can be replaced by 
the distribution function at the equator. The temporal 
evolution of the energy and pitch angle dependence on the ion 
distribution in a 3-D spatial space can be inferred from the 
information on the equator. In solving the kinetic equation 
(1), the time-splitting method is employed. At each fractional 
time step, only one process is considered [Fok et al , 1993, 
1995a]. 

The ring current model has been used to model the recovery 
phases of a moderate storm and a major storm. The model ion 
differential fluxes in general agree well with the observations 
from AMPTE/CCE, except the calculated H + fluxes at energy 
less than 100 keV always exceed those of measurement [Fok et 
al ., 1995a]. However, our model successfully reproduces many 
observed features of storm recovery, for example, drift holes 
on the dayside at energies of 3-20 keV and separation of ions 
as a function of energy due to differences in drift velocities and 
charge exchange lifetimes. The 3-D model also predicts the 
buildup of pitch angle anisotropy observed during storm 
recovery. 



Figure 1, Quiet time equatorial H + fluxes at representative L 
shells [Sheldon and Hamilton , 1993]. 


3. Initial Conditions: Quiet Time Conditions 

We start the simulation with quiet time conditions. The ion 
composition compiled by Sheldon and Hamilton [1993] during 
the quietest days in 1985-1987, seen by the AMPTE/CCE/ 
charge-energy-mass (CHEM) instrument in near-equatorial 
orbit at L- 2-9 R E , is used as initial distribution before storm 
onret. This data set provides average differential ion fluxes in 
an energy range 1-300 keV. The initial energy spectra of 
average H + fluxes at some representative L shells are plotted in 
Figure 1. This quiet time distribution is assumed to be local 
time symmetric, although Sheldon and Hamilton [1993] found 
a small dawn-dusk distortion caused by the cross-tail electric 
field. 

The quiet time PAD of ring current ions is estimated by the 
charge exchange cross sections, assuming that the PAD is 
mainly shaped by the charge exchange loss. The initial value 
of the exponential index n of the pitch angle fit (j(y) = 
j 0 ( 1 +Ay n ), y = sine of the equatorial pitch angle) of H + is 
given as follows, in light of previous modeling results at late 
recovery [Fok et al . , 1995a]: 


« = 73.75 



(3) 


where o s is in 10" 19 m 2 and E is in keV. The value of n 
increases with decreasing L and reaches a value greater than 6 
near the inner edge of the ring current for energies below tens 
of keV. A relatively isotropic distribution of high-energy 
(> 70 keV) H + is obtained, as a result of its large charge 
exchange lifetime. However, using data from the Explorer 45 
satellite, Williams and Lyons [1974] found rounded PADs 
(peak at a = 7t/2) of high energy H + for altitudes in the 
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Figure 2. (top) Dst index and (bottom) the CCE orbit on May 2, 
1986. The times during the two orbit passes in the ring current 
region are 0150-0840 UT for Orbit 1 and 1720-2400 UT for 
Orbit 2. Marks on the orbit are at 1-hour intervals. 

plasmasphere region. They interpreted these PADs as a 
consequence of pitch angle diffusion due to the amplification 
of ion cyclotron waves when the hot ring current plasma 
interacts with the cold plasmaspheric plasma. In section 5, we 
will evaluate the effect of wave-particle interactions on the H + 
PADs, especially at high energies. 

The development during the main phase of a magnetic 
storm, similar to that on May 2, 1986, is studied. The Dst 
index and the CCE orbit on May 2, 1986, are shown in Figure 
2. In the following, we only consider H + , which is the 
dominant ion species for this storm. The simulation starts at 
0200 UT (run time, RT = 0), when the geomagnetic conditions 
are quiet. An empty loss cone distribution is assumed at RT = 
0. During the simulation, an additional loss term with a 
lifetime of half of the bounce period is applied to ions in the 
loss cone, which is defined at 800 km. The measured 
AMPTE/CCE H + flux as a function of pitch angle at L ~ 6.75, 
LT - 4 hours, at RT = 6 (Orbit 1) and 22 (Orbit 2) hours are used 


to update the nightside boundary condition during the storm 
main phase. 

4. Energy Spectra During Storm Main Phase 

The development of ring current H + during the storm main 
phase on May 2, 1986, is simulated with the initial and 
boundary conditions specified in the last section. The offset in 
the convection potential is assumed to be 2 hours, making the 
convection paths symmetric about the meridian of 0800-2000 
LT [Kistler et al. , 1989]. The simulated H + fluxes are compared 
with the CCE measurements during Orbit 2. Figure 3 shows the 
comparison during the inbound (Figures 3a-3c) and outbound 
(Figures 3d-3f) pass at L ~ 5.75, 3.25, and 2.25. There are 
minima in the dayside spectra at about 10 keV seen by CHEM 
measurements. These dips are consistent features of the storm 
time ring current ions, as a result of slow drift velocity 
[Mcllwain, 1972; Fok et al ., 1995a]. As shown in Figure 3, 
our model reproduces the "drift holes" at L < 4. However, the 
simulation predicts a lack of low-energy (below tens of keV) 
ions at the inner edge of the ring current (Figure 3c) due to 
significant charge exchange loss before particles drift from the 
nightside source region to the dayside observation point. 
Significant low-energy H + fluxes of 10 5 s^cm^sr" ’keV* 1 are 
measured. These low-energy ions may be a result of direct 
injection from the ionosphere [Horwitz, 1982] or 
plasmasphere during the main phase. Also radial diffusion 
induced by electric and magnetic fluctuations could diffuse ions 
from the boundary of open and closed drift paths to lower L 
shells. These transport mechanisms are not included in the 
present model. Radial diffusion is believed to be the major 
transport process for high-energy ( E > 280 keV at L - 3) ions 
[Chen et al ., 1993]. However, the good agreement between the 
modeling fluxes and the measured fluxes at high energies 
shown in Figure 3 indicates that this diffusion process has 
only a small effect on redistributing energetic ions in radial 
distance. 

The simulated H + fluxes on the nightside agree well with 
measurements (Figures 3d-3f). At L ~ 5.75, the ion flux is 
similar to a typical geosynchronous distribution, which is a 
double Maxwellian with peaks at a few keV and tens of keV 
[Borovsky et al. , 1994]. At low L shells (Figures 3d and 3e), 
the spectra show the overlapping of two ion populations: the 
freshly injected low-energy ions from the tail and the 
preexistent high-energy particles. In this case, the dips on the 
energy spectra at - 80-100 keV correspond to the transition 
energy between the open and closed drift paths. 

The good agreement shown in Figure 3 between the 
simulated ion fluxes and measurements affirms the ability of 
our model to reproduce energy features of the storm time ring 
current. We next display the spatial distribution of these 
energy features in the average H + fluxes at the equator during 
the main phase of the storm. Fok et al. [1995a] used 
chromograms to display the global energy distribution of ion 
fluxes. A description of chromogram can be found in the work 
by Fok et ai [1995a] and is also given in the appendix. In 
Plate 1, ion fluxes are divided into three energy ranges: 1-5, 
5^40, and 40-300 keV. They are represented by red, green, and 
blue (color wheels in Plate 1) in the chromogram, 
respectively. The geometric mean of the two boundaries for 
each energy range is also labeled on the color wheel. The gray 
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(b) L=3.46-3.03, 1000LT (c) L=2.56-1.97, 1100LT 





Figure 3. The comparison of calculated average H + fluxes (curves) with Active Magnetospheric 
Particle Tracer Explorers (AMPTE)ZCCE measurements (circles) during Orbit 2 at selected 
locations. 


scale bar is a guide to the interpretation of brightness, which 
represents the average flux intensity over the entire energy 
range. Contours of brightness are overlaid on the chromogram 
to distinguish brightness from color. 

Plate 1 shows chromogram representations of H + average 
fluxes at the beginning, middle, and end of the storm main 
phase. The Kp value as a function of universal time on May 2, 
1986, and run time is given in the upper left panel. At RT = 1 
hour (Plate la), the H + flux is mostly a quiet time distribution 
(cf. Figure 1). High-energy H + dominates at L < 3 (dark blue 
area). Low-energy fluxes increase with L shells and are 
indicated by the magenta ring at 3 < L < 5, where high- and 
low-energy ions coexist with a lack of medium-energy 
particles. An injection of ions is seen at the nightside 
geosynchronous region. The injections are mainly carried by 
ions of energies less than 40 keV (yellow fringe). In response 
to the intensifying convection field, the freshly injected ions 
move earthward and result in a strong day-night energy 
asymmetry in flux level at RT = 9 hours (Plate lb). The total 
ion flux is increasing during the main phase. It should be noted 
that the intensity range indicated by the gray scale bars 
increases for each successive panel in Plate 1. As shown in 
Plate lb, less energetic ions corotate and reach the day side 
through dawn (red area in the dawn-noon quadrant). Energetic 
ions drift westward and may move out of the model boundary in 
the noon-dusk quadrant (green-cyan area). The blue area near 
the inner edge of ring current is the preexisting high-energy 
ions, whose drift paths are closed so they cannot be 
transported from the tail through open drift paths. The green 
and blue areas form a boundary between open and closed drift 
paths for the low-energy (< 5 keV) ions. At RT = 21 hours 


(Plate lc), Kp reaches a value of 7 + . Ions penetrate closely to 
the Earth on the nightside. The maximum ion flux is located at 
L ~ 4.5. The light yellow represents a flat energy distribution 
of ions in this region, with a slight peak at low energy. The 
area forbidden to the convecting low-energy ions is 
compressed in response to the strong convection field. 

5. Pitch Angle Distribution 

The spin of the CCE spacecraft allows particle sampling at 
all pitch angles [Williams and Sugiura , 1985]. Figure 4 shows 
the measured differential fluxes (circles) in four energy ranges 
as a function of the equatorial pitch angle during Orbit 2. Data 
in the inbound (dayside) pass are shown except for 
measurements in the energy range of 17.9-28.2 keV, where 
data in the outbound (nightside) pass are plotted because of the 
low flux level on the dayside at these energies. At L ~ 2.75 
(left panels), ion fluxes of all energies strongly peak at 90° as 
a result of rapid charge exchange loss of particles with small 
pitch angles. The distributions are flatter as L is increasing, 
where the neutral hydrogen density is low. Near 
geosynchronous orbit (right panels), peak fluxes of low- 
energy (< 30 keV) ions are found at field-aligned pitch angles. 
This enhancement of field-aligned low-energy particles during 
storm injections has been observed by the ATS 6 
geosynchronous satellite and has been suggested as a signature 
of the parallel energization mechanism [ Mauk , 1986]. 

Simulated PADs plotted in solid curves are overlaid with the 
CCE measurements in Figure 4. Calculated fluxes are scaled by 
the measured fluxes near 90° pitch angle in order to compare 
the shapes of the distributions. The calculated field-aligned 
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Plate 1. Chromogram representation of the average flux of H + in s^cnrf-^sr^keV" 1 at RT equal 
to (a) 1, (b) 9, and (c) 21 hours. 
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pitch angle (degree) 


Figure 4. The comparison of calculated H + fluxes as a function of the equatorial pitch angle 
(curves) with AMPTE/CCE measurements (circles). Dashed curves in the bottom panels (100-154 
keV) represent distributions when an additional pitch angle diffusion, with D aa = 5xl0“ 6 s' 1 , is 
considered. 


distribution at 1.47-2.19 keV, at L ~ 5.75 (top right panel) 
results from the boundary PAD at L - 6.75. It can be seen that 
our model agrees with data very well at low energies. This 
implies that PADs of low energy ions are mainly controlled by 
the drift motion and the charge exchange process, which are 
considered in our model. At high energies (bottom panels), the 
model predicts a flat distribution of trapped ions and a sharp 
drop-off near the loss cone. This predicted flat distribution on 
the dayside is a result of the isotropic boundary condition on 
the nightside. However, the data on the dayside show a gradual 
decrease in flux from the trap region to the loss cone. It seems 
that some diffusion processes in pitch angle, most likely due 
to wave-particle interactions, must take place to smooth out 
the distributions. Kennel and Petschek [1966] showed that for 
ion energies greater than E c = B 2 l2ji Q n b , where n b is the cold 
plasma density, the interaction with cyclotron waves leads to 
pitch angle diffusion. At L = 5, assuming n b = 10 cm' 3 , E c ~ 15 
keV is obtained. This estimate supports the argument that 
interactions with ion cyclotron waves play a more important 
role in determining PAD for ions of energies at or above tens 
of keV than ions of lower energies. 

In order to estimate the wave amplitude that could account 
for the measured PAD, an additional term of pitch angle 
diffusion is added in the kinetic equation (1), with the diffusion 
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Therefore { D KK ) can be calculated for a given D aa . Since the 
phase space density in our model (1) is expressed in terms of 
(M y K)y at each time step, before the diffusion term of (7) is 
solved, f s at each K value is mapped into a fixed E grid. After 
the diffusion in K is performed, f s is mapped back into M 
grids by interpolation and then other processes, such as drifts, 
charge exchange, and Coulomb losses are performed. The 
resulting PADs with this additional diffusion process for 100- 
154 keV ions are shown by dashed curves in the bottom panels 
of Figure 4. In this case, D aa = 5xl0" 6 s' 1 provided a good 
match to the data. 

Next we estimate the wave intensity corresponding to the 
diffusion coefficient obtained above. Assuming a pure pitch 
angle diffusion, the wave frequency is much less than the 
gyrofrequency of the ring current ions. In the cold plasma 
theory, the diffusion coefficient can be roughly approximated 
as [Cornwall et ai , 1970] 


®aa 


SB 

B 


2 


(9) 


where (o g is the ion gyrofrequency and SB is the total wave 
amplitude. For 100 keV H + at L = 5, a D aa of 5xl0" 6 s' 1 gives 
SB -0.1 nT. This amplitude of magnetic field fluctuation is 
consistent with the low-frequency wave observations by the 
AMPTE [LaBelle et aL , 1988; Anderson et al., 1992] and the 
GEOS satellites [Roux et al.y 1982] and corresponds to a 
moderate ion cyclotron fluctuation. 


6. Induced Electric Field 

It is possible to obtain a qualitatively realistic ring current 
from a simple global enhancement of the magnetospheric 
convection electric field [Chen et al.y 1994]. However, it is 
well known that the nightside magnetic field undergoes a 
series of stretching and relaxation cycles, each associated with 
isolated substorm events [e.g., Moore et al. , 1981], that each 
of these cycles is associated with a “convection surge” [Quinn 
and Southwoody 1982], and that particle acceleration 
associated with these surges contributes in some way to 
magnetic storm periods [Mauk and Meng y 1983]. The induced 
electric field associated with these events is not global but 
rather localized in the evening sector. Lyons and Williams 
[1980] argued that this localization of the injection fields is 
quantitatively critical to the formation of the ring current. 
Chen et al. [1994] constructed the substorm-associated 
convection surges by a superposition of randomly occurring 
impulses that rise sharply and decay exponentially with a 
lifetime of 20 min, but treated these as variations of the cross- 
tail potential drop or, equivalently, the global convection 
strength. 

In an attempt to incorporate the localized character of 
substorm-induced electric fields, we have adopted an approach 
that is grounded in observations of plasma dynamics at 
geosynchronous orbit, which corresponds well with the outer 
boundary of our simulation space. To a steady global 
convection pattern, we add the drift velocity (ExB) induced by 
changes in a model magnetic field, using an approach that is 
similar to that used by Delcourt et ai [1990] to describe 


variations of the fields in the more distant tail. That is, the 
Tsyganenko magnetic field model (in this case, Tsyganenko 
[1989] EXT89AE routine) is given a time dependence through 
the observed variations of the AE index, plotted at 1-min time 
resolution in Figure 5 for the day of May 2, 1986. The 
magnetic field configuration is associated with six levels of 
the AE index, as shown in Table 1. Our procedure honors this 
association of the Tsyganenko level with AE, while 
introducing large variations around the average at moderate 
activity levels, as follows. 

The Tsyganenko level (TL) is initially set at a level, 7X,, 
that could be determined by knowledge of the instantaneous 
geotail configuration, if known, but is treated as a free 
parameter. TL is assumed to closely track AE at low values of 
AE, according to the Tsyganenko relationship. However, it is 
taken to grow linearly at a rate proportional to the excess of 
AE over a critical value, AE 0 . The growth rate tracks AE in such 
a way that the induced flow increases with AE. To simulate 
collapse of the tail owing to an instability, TL abruptly 
decreases immediately upon reaching a threshold value of 
7X max (generally set equal to 6, the maximum TL available). 
The amplitude of the decline depends on the AE value at the 
beginning of the field line collapse. However, the minimum 
duration of the expansion phase is set to be 2 min. After each 
collapse or expansion phase, TL then resumes tracking the AE 
index as before. During the growth phase, earthward 
convection in the tail is suppressed by this procedure, while 
during the expansion phase it is dramatically increased. The 
average rate of convection increases with AE at low levels, 
while tracking it closely. At higher levels, convection 
undergoes cycles of suppression and enhancement, but the 
average convection rate must ultimately match that of the 
underlying convection, a basic constraint that this procedure 
preserves. 

This algorithm for modeling the induced electric field as a 
series of substorms is summarized mathematically in the 
following: 

Initial level (free) 

TL = TL i (10a) 

Growth or recovery 

^ = C l (AE-AE 0 ) (10b) 


Expansion 



duration = max 


2 min, 


2(6- TL e ) 

Cl 


(10c) 


where AE 0 = 25 y, C\ - 2x1 0" 6 s" 1 ^ 1 , ^2 = 0-25 min" 1 , and TL e 
is the Tsyganenko level according to Table 1 at the beginning 
of the expansion. C\ is chosen in a way that the ion drift on 
the nightside in response to the field line stretching will not 
exceed the earthward convective drift to prevent net tailward 
motion on the nightside. C i is chosen such that the complete 
dipolarization from level 6 to 1 takes place in 20 min. The 
change in Tsyganenko level during dipolarization is 
controlled by the duration in such a way that TL e is at the 
middle of EL max and the bottom level at the end of the 
dipolarization. However, during very active periods that 
whenever TL e > 5.75, the duration of dipolarization is set to be 
2 min and TL drops from level 6 to 5.5. 

Following the above algorithm and with TL t set to be 5.5, 
the simulated Tsyganenko level as a function of time during 
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Figure 5. (top) 1-min AE index on May 2, 1986. (bottom) Simulated instantaneous magnetic 
configuration represented by levels of the Tsyganenko model. 


the main phase of the storm is plotted in Figure 5 together 
with the AE index. As shown in Figure 5, there are five distinct 
large-amplitude dipolarization events at about 0320, 0730, 
0855, 1025, and 1220 UT. After 1220 UT, the geomagnetic 
activity is so strong that repeated rapid stretching and 
collapsing of the magnetic field occurs accompanied by only 
small-amplitude changes in TL. The predicted substorms 
qualitatively resemble the injections of energetic electrons 
observed by the Los Alamos National Laboratory (LANL) 
geosynchronous satellites. LANL data for May 2, 1986, show 
clear particle injections at about 0300, 0700, 1030, and 1220 
UT, and injection cycles are nearly continuous after 1220 UT 
(G. Reeves, unpublished data, 1995). Although the model does 
not exactly reproduce the number or timing of injections 
reported by the LANL satellites, it duplicates the important 
features of isolated substorms before 1220 UT and multiple 
substorms in rapid succession thereafter. The free parameter 
TLi was chosen to obtain good agreement in time between the 
predicted and the observed discrete injections. 

The most direct approach to calculating the drift of ions as a 
result of the changing magnetic field might be to calculate the 
induced electric field by “uncurling’* <7B Idt and then calculate 
the corresponding EXB drift. However, when using an 
empirical magnetic field model, it is considerably more 
straightforward to apply Alfven’s theorem, also known as the 
frozen-in field theorem. For purposes of computing the induced 
field motion, each field line can be regarded as fixed and rooted 
in the ionosphere, since magnetic field variations are 


vanishingly small there. The EXB drift velocity of an ion at 
any point along a field line is equivalent to the perpendicular 
velocity of that field line at that particular position [Wolf, 
1983]. Therefore the induced equatorial drift velocity is 
calculated from the instantaneous velocity of the equatorial 
crossing of magnetic field lines in response to the temporal 
variation of magnetic level shown in Figure 5. This substorm- 
induced drift is added to the gradient curvature drift and the 
convection electric drift to move the ions. Though the 
motions of the Tsyganenko magnetic field lines during 
substorms are considered in calculating the induced drifts, a 
constant dipole field is assumed in the rest of the calculations. 
This is a rough and non-self-consistent way of implementing 
the time-varying magnetic field, which is adequate only in the 
inner magnetosphere. A more self-consistent treatment will be 
pursued in future work. 


Table 1 . AE parameterization of 
Tsyganenko Magnetic Field Level 


Level 

AE 

1 

0-50 

2 

50 - 100 

3 

100 - 150 

4 

150 - 250 

5 

250 - 400 

6 

400 and up 
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(a) quiescent radial convection at Kp = 7 


(b) induced radial drift during dipolarization 
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Figure 6. Radial components of (a) steady convection drift velocity at Kp = 7 and (b) induced 
velocity during dipolarization. 


Figure 6 compares the radial component of quiescent 
convective drift at Kp = 7 with that of the induced velocity 
during dipolarization. The corresponding azimuthal electric 
field that cause these drifts can be estimated by E<p = 5- v rad ia] . In 
this case, the induced electric field at midnight at L ~ 6 during 
dipolarization is about 1.3 mV/m, which is comparable with 
the quiescent convection field at the same location. However, 
the inductive electric field is more localized, having E ^ « L k , 
where k ~ 1.2, compared with k = 1 in the quiescent field 
model. As shown in Figure 6, the magnitude of the steady 
convection is day-night symmetric (Figure 6a). Ions with open 
drift paths moving in from the tail will drift out the dayside 
magnetopause with the same speed. Most of the energy gained 
by particle injection and energization on the nightside is lost 
by particle drift out at the dayside magnetopause [Kozyra et al, 
1994]. However, during the main phase of a storm, the steady 
intensification of convection leads to deeper penetration of 
ions (Plate la) and increase of particle energy. With the 
consideration of changing magnetic field, the drift due to the 
induced electric field is stronger on the nightside than on the 
dayside during the expansion of a substorm (Figure 6b). The 
large resultant velocity of convection and dipolarization on 
the nightside pushes ions close to the Earth before they are 
trapped [Lyons and Williams , 1980]. As a result, the energy 


content in the ring current has a sudden jump during the 
dipolarization event. The situation during the field line 
stretching is the opposite. The induced radial drifts during 
stretching are in the opposite directions of those shown in 
Figure 6b, and with smaller magnitudes. The tailward motion 
of field lines counter balances the steady earthward convection 
on the nightside. The drift loss at the dayside magnetopause 
thus exceeds the input from convective energization. The 
energy content in the ring current may drop during the 
substorm growth phase. 

Figure 7 shows the H + energy as a function of simulation 
time with and without the inclusion of the induced electric 
field. As shown in the figure, the quiescent convection causes a 
steady growth in the H + energy (solid curve). In the case of an 
additional inductive drift, the H + energy fluctuates as field lines 
undergo stretching and collapsing. During the active period of 
10 to 22 RT, the H + energy oscillates rapidly about a mean 
energy that is a bit less than the energy obtained by steady 
convection only. It was initially very surprising to us that less 
ring current is produced when storage-release cycles of 
localized inductive convection are superposed upon a steady 
convective flow. Upon closer examination of the results, we 
see that this is true only for the early main phase, when 
convection is weak so that large cycles of growth and 



Figure 7. H + energy as a function of simulation time during the storm main phase with and 
without the consideration of the induced electric field. 
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dipolarization occur. During that phase, inclusion of the 
inductive fields clearly fails to produce as much ring current as 
the steady case with the same mean convection. However, 
during the later phase, when convection has increased so that 
growth is as fast as dipolarization, we see that the two cases in 
Figure 7 track each other closely. 

One interpretation of this is that localized sporadic 
injection is not an efficient process for creating a full ring 
current. At low storm power levels, in the absence of a strong 
steady convection, the freshly injected particles are not able to 
form a symmetric ring current. As a result, the bursty increase 
of the ring current during dipolarization is more than balanced 
by the reduced convection during the slow and long (on the 
order of hours) growth phase. At higher storm power levels, 
the cycles of growth and dipolarization are small-amplitude 
perturbations on the fast global convection that has become 
the dominant process. The high power levels smooth out the 
response of the magnetotail, which spends most of its time in 
more highly stretched states, and transmits energy into the 
ring current relatively directly and without significant storage 
and release. 

The large amplitudes of the substorm-associated fluctuations 
in the ion energy content shown in Figure 7 are not supported 
by the H component of the high time resolution (1 min) 
symmetric disturbance index (SYM-H) [ Iyemori , 1990], which, 
on May 2, 1986, fluctuated with amplitudes smaller than those 
calculated (T. Iyemori, unpublished data, 1995). However, 
direct inference of ring current energy content from the surface 
field perturbation requires steady configurations of the 
magnetic field and the trapped particles [ Sckopke , 1966]. A 
more accurate way to predict the magnetic perturbation is to 
calculate the current density produced by the energetic ions. It 
was found that the pressure gradient is the main contribution to 
the current density [Lui et a/., 1987], and the effect is scaled 
inversely by B. During substorm injections, particles are 
transported earthward and experience a stronger magnetic field. 
The increase in the current density produced by these injected 
particles may somehow be reduced by the larger local B values. 
We expect the field depression produced by the ring current to 
fluctuate with a smaller amplitude than the energy content does 
during substorms when the self-consistently time-dependent 
magnetic field is taken into account in the calculation. 

7. Discussion and Summary 

The evolution of the ring current ions from an initially quiet 
time distribution has been addressed as a response to the storm 
time enhanced convection field and particle injection. 
Although simple models for the magnetic and electric field are 
used, our model reproduces the major features of the storm time 
ring current, such as the day-night asymmetry of ion 
distributions and deep penetration of ring current ions at the 
peak of the storm. Under the assumption of a magnetic dipole 
field, there are simple analytical expressions for the bounce- 
averaged drift velocities [Ejiri, 1978]. Others have calculated 
the drift motion in a nondipole field [Schulz and Chen , 1995; 
Chan et al , 1995]. 

We have introduced a scheme for solving pitch angle 
diffusion in an M-K space. Phase space densities are expressed 
in terms of adiabatic variables M and K in solving drifts and 
collisional losses, while E-K parameterization is used for the 
pitch angle diffusion. If an M-K parameterization is used 
throughout the calculation, a pure pitch angle diffusion may 


result in diffusion in M , diffusion in K, and terms 
corresponding to mixed diffusion of M and K. The advantage of 
switching between two domains is that one can have a simple 
expression for the left-hand side of the kinetic equation (1), 
without terms of M and AT , and can simultaneously avoid 
dealing with mixed diffusion on the right-hand side of the 
equation. 

The effect of the stretching and collapsing of the 
magnetosphere is represented by an additional drift velocity 
induced by this field line motion. As a consequence, the energy 
content in the ring current oscillates about a mean value, 
which does not change significantly because of the substorm- 
induced activities but grows steadily as the convection electric 
field is enhanced. This result supports the two-step process of 
ring current formation: during dipolarization events, the 
earthward transport of plasma sheet population forms a 
spatially sharp boundary of energetic plasma near the 
geosynchronous orbit [Mauk and Meng , 1987], then particles 
undergo further earthward displacement and are energized 
adiabatically by the convection electric field. The time scale of 
substorm activities is on the order of 10 min or less, which is 
shorter than or comparable to the bounce periods of low- 
energy ions (< 1 keV) at L > 5. Under this condition, the 
assumption of conserved J (or K) and the concept of bounce 
averaging are no longer valid. The lower limit of the particle 
energy has to be high enough to ensure the conservation of the 
second adiabatic invariant. 

In summary, we have used a 3-D ring current model to 
simulate the main phase of a magnetic storm similar to that on 
May 2, 1986. Modeling results are compared with the 
AMPTE/CCE particles measurements. We found 

1. The calculated average H + fluxes agree well with 
measurements, except the model cannot reproduce the observed 
ion fluxes at energies below tens of keV at L < 2.5. 

2. Near the inner edge of the ring current, ion flux peaks at 
90° pitch angle, as a result of strong charge exchange loss of 
the field-aligned ions. 

3. For H + energies greater than tens of keV, the round shape 
of the pitch angle distribution cannot be solely explained by 
the charge exchange process and the drift motion. Diffusion in 
pitch angle due to interactions with ion cyclotron wave are 
suggested to be responsible for the observed distribution. A 
pitch angle diffusion process with diffusion coefficient about 
5x1 O' 6 s' 1 is required in order to match the data and is 
consistent with observed fluctuation amplitudes. 

4. The energy content in the ring current fluctuates in 
response to the field line stretching and collapsing during 



Figure Al. Curve j(E) is divided into three subranges that are 
associated with red, green, and blue light. 
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Table A I. Guide to the Interpretation of Chromogram Color 


Color 

Interpretation 

Gray 

no dependence on E (flat) 

Red 

strongly peaked at low values of E 

Green 

strongly peaked at midscale in E 

Blue 

strongly peaked at high values of E 

Yellow 

weakly peaked toward low E 

Cyan 

weakly peaked toward high E 

Magenta 

peaked at high and low E 


substorms. However, the overall increase in the ring current 
energy is a result of the gradual enhanced convection field. 

Appendix: Description of a Chromogram 

The goal of a chromogram is to plot a function of three 
variables, j(L, <f>, E), in a two-dimensional color plot. The 
dependence on the First two variables (i.e., L and <p) forms a 
pixel image in a two-dimensional L-<p space, while the 
dependence on the third variable (E) is encoded as color. The 
encoding of color uses the concept of three basic colors 
(photon energies) of light; red, green, and blue. Variable E is 
divided into three ranges (Figure Al), associated with red (R), 
green (G), and blue (B), respectively. The average values of j in 
these ranges are used to scale the R, G, B components of the 
encoded color. The hue of the resulting color corresponds to 
the centroid or mean E of j(E). The saturation of the resulting 
color corresponds inversely to the width of j(E). Saturation 
measures the degree of dominance of one hue over the others. 
As color is varied from gray to a pure hue (from the center to 
the perimeter of a color wheel shown in Plate 1), saturation 
varies from 0 to 100% at constant brightness. The brightness 
of the encoded color corresponds to the average j over the full 
range of E . Table Al summarizes the interpretation of the color 
code. The rough energy distribution at each L-0 location can be 
inferred from the color at that particular pixel. The detailed 
features of energy spectra cannot be resolved using the 
chromogram representation, but it is useful to indicate the 
spatial dependence of the energy distribution. 
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